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Mechanistic Model to Replace Hodgkin-Huxley Equations 

Bo Den^ 

Abstract: In this paper we construct a mathematical model for excitable membranes by in¬ 
troducing circuit characteristics for ion pump, ion current activation, and voltage-gating. 
The model is capable of reestablishing the Nernst resting potentials, all-or-nothing action po¬ 
tentials, absolute refraction, anode break excitation, and spike bursts. We propose to replace 
the Hodgkin-Huxley model by our model as the basis template for neurons and excitable 
membranes. 

1. Introduction. When Hodgkin and Huxley eonstrueted their model for the squid giant axon ('ifTIl 
they were fully aware of their model’s drawbaeks beeause it was only a phenomenologieal fit to 
their experimental data. They eommented speeifieally that different empirieal forms should fit the 
same data or even better. They were right and other researehers saw the same problem too [|2l[3l|4l. 
Alternative models were eoneoeted (dll) but never gained any traetion beeause there was no point 
to replaee one ad hoc model by another arbitrary one. However, replaeing a phenomenologieal 
model by a meehanistie one is a different matter entirely. 

By meehanistie it is meant for a model to have as few hypotheses as possible that apply to phys- 
ieal proeesses or objeets of a same type. Newton’s inverse-distanee-squared law for gravitation is 
the first and one perfeet example of meehanistie modeling beeause it applies to all maeroseopie 
bodies of mass. Goldman’s derivation (lEl) of Nernst potentials aeross eell membranes is meeha¬ 
nistie beeause it applies to all ion speeies. In eontrast Hodgkin and Huxley’s individual treatments 
of the sodium and the potassium eurrents are not meehanistie beeause their hypothesis for the 
sodium ion does not apply to the potassium ion or viee versa. Researehers must have tried but 
failed beeause other than various variations of the HH model no meehanistie model ean be found 
in the literature. 

The purpose of this paper is to fill this literature gap. The idea is to model the membrane 
as a eireuit of deviees eaeh is defined by a eurrent-voltage eharaeteristies. We will model the 
sodium-potassium ion exehanger pump by the IV -eharaeteristies that the time-rate of ehange of the 
eurrent is proportional to the power of the pump. We will model the ion ehannel aetivation and the 
voltage-gating by one unified IV -eharaeteristies that the voltage-rate of ehange of the eonduetanee 
is proportional to the eonduetanee. We will demonstrate that the resulting eonduetanee-adaptation 
model is eapable of reprodueing all known phenomena of the HH model and mueh more, and 
henee provides a meehanistie alternative to the HH model and a model template for other types of 
exeitable membranes in general. 
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2. The Result. The approach to modeling the squid giant axon as an electrical circuit began 
with Cole’s work (d^ 111 IH) that the axon membrane had a rather narrow range around l/iF per 
square centimeter for membrane capacitance. Let V be the intracellular potential difference (inside 
potential minus outside potential) and I = CdV/dt denote the capacitance current. Let Jk, Ins. 
denote respectively the potassium and the sodium ion current using the outward direction as the 
reference direction. Rather than the leakage current considered by the HH model we consider 
instead the voltage-gating current Iq measured outward. Last let Jext be the net remaining current 
with the inward reference direction. Then the following equation must be the first equation of any 
circuit model of the membrane by Kirchhoff’s Current Law 

dV 

c = — [Jk + ^Na + -^G “ -^ext]- 

The aim of this section is to derive the functional forms for the potassium, the sodium, and the 
gating currents based on our proposed mathematical models for ion pump, ion channel activation, 
and voltage-gating, and to obtain the resulting model for the membrane. 


Ion Pump. Let q denote the intracellular charge difference (charge concentration inside minus 
outside) of a given ion species across the cell membrane. The following as a mathematical model 
for the pump was proposed in @ 


dt'^ 


, dq 

V 


( 1 ) 


with 0 a parameter. If we let / = dq/dt he the current through the ion pump and V = q/C he 
the across-membrane potential over the pump with C being the membrane capacitance. Then the 
model can be interpreted as the /^-characteristics of the pump as an electrical device because the 
model is equivalent to 


dt 


XIV 


with A = C(/) referred to as the pump parameter in the unit of per time per voltage. As the product 
IV represents power. The pump characteristics can be stated as the change of the ion pump current 
is proportional to the power of the pump. This is a way to model the energy transfer of the pump 
when ATP is converted to ADP in exchange for ion transportation across the membrane. 

The minimal circuit consisting a pump, a resistor, and a capacitor in loop can be understood 
completely, c.f. Fig{TJ The ordinary differential equations for the circuit is 


CV' = I, J'=-A/(1/ + 7 /). 


( 2 ) 


Eliminating the time variable the equation becomes dl/dV = —XCiV + 7/) which being a linear 
equation can be solved explicitly. In fact all important insights can be obtained by dropping the 
resister term (7 = 0), for which the solution with the initial condition V (0) = 0, /(O) = Jq > 0 is 


V(t) 


1 _ p AV^t p AV^i 

Y ^ — o\/^T/2 ^ 


l + e 


-XVoot 


and I{t) = 2XCV^ 


[1 -f e->'Voot'j2 


( 3 ) 
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Figure 1: (a) Circuit diagram of the potassium ion pump with the capacitive membrane, (b) The 
phase plane portrait of the circuit equations dH) with 7 = 0 . (c) A CRE-cixcmi approximation 
of the pump circuit, (d) The Nernst potential under electrical neutrality and across-membrane 
diffusion. 


with 



The solution for the current means the pump is always unidirectional. That is, if it is positive 
at one point, say Jq > 0, then it stays positive all the time. Equally important, as f —)■ —00, 
liml/(f) = —Vqo, lim/(f) = 0 and as f ^ +00, liml/(f) = 14o, lim J(f) = 0. This means if the 
membrane is negatively charged in the past (~ — 14o < 0) then the membrane will be positively 
charged eventually (~ 14 o > 0) if the pump’s transporting direction is from outside to inside as 
with the case of the potassium pump. Similarly, for the sodium pump, it will transport all sodium 
ion from inside to outside in asymptote. 

We demonstrate next that the minimal pump-capacitor circuit is exponentially close to a resistor- 
battery-capacitor circuit in loop with the battery value given hy E = — 14o, c.f. Fig Ill The differ¬ 
ential equation for the circuit is: CV = I = —{V + E)/r with the resistance r to be determined. 
Solving it with the zero voltage initial condition 1/(0) = 0 as we did for the pump circuit we have 


V{t) = -E[l - = l/oo[l - 
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So if we choose the resistance as 


( 4 ) 


then the dynamics for the battery circuit and the pump circuit are exponentially close for large t 
since the pump dynamics is approximately l/(f) = 14o[l—~ 14o[l— 
converging to the same resting potential Voo at the same exponential rate and starting from 

the same initial value V (0) = 0. That is, the ion pump can be approximated by a conductor-battery 
in series with the conductance 1/r = AC14o. This approximation becomes more accurate as time 
goes by. This gives a mechanistic justification for Hodgkin and Huxley’s modeling of ion channel 
currents by battery-driven ion channels. 


Nernst Potential. The pump dynamics above says that in the absence of other forces such as 
diffusion a particular ion pump will eventually deposit all ions of the same type from one side of 
the membrane to the other side. We determine next the value of the asymptotic potential 14o of 
the pump under the influence of diffusion and the electrical neutrality condition of the membrane, 
c.f. FiglU which is an illustration for the potassium ion. By Bernstein’s hypothesis (' [[TOll l this 
inward current due to the electromotive force is balanced out by the outward diffusion assuming 
the membrane is permeable to potassium ions. Specifically, let Ip, denote the pump current 
and the diffusive current as shown. Let the reference direction be inward and be defined by the 
across-membrane variable x from outside to inside and a be the thickness of the membrane. Let 
n{x) be the ion density at position x and let Ek = —Voo be the equilibrium state reached when 
these two currents cancel each other out under the neutrality assumption L = 0. Then the same 
Bernstein-Goldman equation below must hold ([[SinTIl 


0 — Id + Ip 






with D the diffusion constant and /i the mobility parameter (' [[TTll l. Solving this linear equation 
in n from a; = 0 to a: = a, and then expressing the result in Ek to obtain the following which is 
exactly the potassium ion’s Nernst resting potential 


Ek = 



]K]i 

1K|. 


with [K]j = n(a), [K]o = n(0) denoting the inside and outside potassium concentration, respec¬ 
tively. The same derivation can be used to obtain the sodium ion’s Nernst potential. 


Ion Channel Activation. With the presence of various types of ions and charged molecules, the 
across-membrane potential V is the manifestation of their aggregate. The Independence Hypoth¬ 
esis in electrophysiology holds that the opening and closing of an ion gate is a function of the 
voltage not of the other ions. As a consequence individual ion currents are modeled by the ohmic 
/^-characteristics form I = g(V — E). By eliminating the spatial effect of axon by their voltage- 
clamp experiments Hodgkin and Huxley showed l lfT^ [T3l [T4l [TSll l the conductance g undergo 


4 




changes throughout the course of an aetion potential and had the insight to propose it to be a fune- 
tion of the voltage, that is, = g{V). Their work showed that throughout the course of an action 
potential the ion eonduetanees for sodium and potassium ions, inerease from near zero 

when the membrane is at rest and then deerease to near zero again when the aetion potential ends, 
henee implying that there are ion ehannels whieh open to inerease their eonduetanees and elose to 
deerease them. Here is where the HH model beeomes empirieal rather than meehanistie beeause 
of its use of arbitrary functional forms for the sodium and potassium ehannel eonduetanees. 

By ehannel activation we mean channel opening and elosing. We propose the following model 
for both types of aetivation 


^g g^V equivalently 


dg g 

dV b 


( 5 ) 


with 6 > 0 the aetivation parameter in the dimensional unit of voltage. This means V is inereasing 
(Al/ > 0) if and only if the conductanee g is inereasing (Ag > 0). That is, the ehanging of 
eonduetance (opening or elosing of ion channels) with respect to depolarizing or hyperpolarizing 
potential is proportional to the eonduetance. It is a double-edged sword — positive feedback with 
increasing voltage and negative feedbaek with the opposite. Solving it we get the exponential 
aetivation law g{V) = . Normalizing it at the ion’s Nernst potential with g[E') = g we get 




Namely, a = ge~^^^. The parameter g{E) = g is referred to as the intrinsie eonduetanee. It is 
a simple exereise to eheek that the eorresponding IV-cmve I = fiV) = ge^^~^db(y _ 
a unique minimal point atb = E — b whose value is —gbe~^, e.f. Figj2l To the right side of b, 
the /(/-curve / is increasing, erossing the V axis at the only resting potential E. To the left side 
of the minimum eriticality b, f is deereasing, giving rise to a varying negative conductanee. It 
approaehes 1 = 0 asymptotieally from below as V approaehes negative infinity. For reasons whieh 
will become clearer later we will refer to parameter b as the activation range parameter, measuring 
the minimal eurrent point from the Nernst potential. 

The proposed aetivation eharaeteristies has this property of diehotomy. If the membrane volt¬ 
age is inereasing in time, dV/dt > 0, then the aetivation gates open up exponentially fast in 
eonduetanee as dg/dt = dg/dV ■ dV/dt = kg with k = (dV/dt) // > 0. In eontrast, if the mem¬ 
brane voltage is deereasing in time, dV/dt < 0, then the aetivation gates elose down exponentially 
fast as well as dg/dt = kg with k = {dV/dt)/b <0. It is a model for both opening activation 
and closing activation with the same voltage-speeific exponential rate. That is, opening gates beget 
more gates opened and closing gates triggers more gates elosed. 

For the potassium and sodium /(/-characteristics we get respeetively 

/k(V') = {V - Ek), f„.(V) = g„.e^'^-‘^-‘>'‘HV - Bn.). 
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(a) (b) 


Figure 2: (a) Typical Jl/-characteristics curves for the sodium, potassium, respectively. 

The left-most intersection of I = and / = is a first approximation of the 

membrane resting potential Em- But the geometric configuration of the two Jl/-curves implies 
the KNa-equilibrium state i^KNa is almost always unstable. Here E'Na — is the minimum point 
of the -curve, (b) A similar plot as (a) except for the addition of the gating IV-curve with 

Eq = £^KNa- The true resting membrane potential Em is the left-most intersection of J = 
and I = —(/n^(I^) + shown in the inset. 


In the operating range of action potentials between the potassium and sodium Nemst potentials, 
{Ek, T^Na), only the critical potential *^he sodium current may lie, predicting that the 

potassium current is always outward and the sodium current is always inward which in turn has two 
different phases: increased entry into the cell to the left of the critical point i?Na —^Na decreased 
entry to the right of the criticality. That is, inside this operating range the sodium and potassium 
IV-curves behave differently qualitatively, the former usually have a negative conductance branch 
but the latter’s conductance is always positive. 

There is a quantitative difference between the intrinsic conductances and First we know 
from the derivation of the Nernst potential above (specifically (@1)) that the conductance is 

9,(0) = 9,6-®""’'' = i 

with Voo ~ |.E'k|. There is no need for a precise value of 14o, an estimate of its range or average of 
its range suffices. Say E^ ~ —60 and V^o ~ 60. This gives rise to an estimate of 

g^ ~ AC60e^K/6K = AC60e-®°/^K. 

Similarly, if we use £'Na ~ +70, a similar estimate is obtained for the sodium conductance: 

g^^ ~ ACTOe-^N.+N, = AC70e^°/'’Na. 
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The conclusion is if the 6-parameter values for both ions are comparable, say in the decade range, 
then relative to each other, the intrinsic potassium conductance should be a few orders smaller 
than the intrinsic sodium conductance From the HH model the intrinsic potassium conduc¬ 
tance is worked out to be g^ = 0.0229 m.mho per square centimeter. As a result a guessed value 
is used for the intrinsic sodium conductance g^^ = 100 m.mho/cm^, both are used below and for 
the illustrations of the IV -curves in Figl2l 


Voltage-Gating. A major advance in neurophysiology was made in the 1970s by the discov¬ 
ery of the voltage-gating phenomenon whereby there is a small pulse-like outward current opposite 
to the inward entry of sodium ions during the onset of action potential across the squid giant axon 
membrane when the membrane is depolarizing. The gating current is due to the release of charged 
molecules from the sodium channel pores in responding to some conformational changes of the 
pores to the depolarizing voltage. (This gating current was not recognized by Hodgkin and Huxley 
but was nonetheless captured by their meticulously fitted model.) 

For the IV -curve of the gating current we propose first it has a resting potential Eq, which may 
be taken to be the intersection of the sodium and potassium /(/-curves, calling it the provisional 
resting potential i^KNa. and second, it is nonlinear ohmic of a similar form as for channel activation 
I = /g(1/) = gG{V){V — Eq) except for a negative proportionality: 


dg g 

dV ~ b 


( 6 ) 


with 6 > 0. That is, positive proportionality is for activation and negative proportionality is for 
gating. Gating conductance recedes in proportion to itself with depolarizing voltage (Ag < 0 if 
and only if AV > 0). Solving it the corresponding /(/-curve becomes 


I = MV) = - Ea). 


For this type of function, it has the unique maximum point at the critical voltage Eq -f right 
of the gating equilibrium Eq. For Eq < V < Eq -f b^ the gating current is increasing in V 
and outward, preventing the membrane from depolarization. Its maximum defines the minimum 
threshold for excitation currents to clear in order for action potentials to generate. For this reason, 
parameter b^ is referred to as the gating range parameter, similar to the activation range parameters 

K.- 

Notice that once the threshold is cleared and the membrane is depolarizing (below V = 0 and 
dV/dt > 0), the gating conductance drops exponentially fast with dg/dt = dg/dV ■ dV/dt = —kg 
with k = {dV/dt)/b > 0, permitting the generation of an action potential. When the membrane is 
hyperpolarizing (below V = 0 and dV/dt < 0), gating becomes active again with the conductance 
approaching the intrinsic conductance g^. To the left side of the gating reversal potential, V < Eq, 
the gating current I = ((/) is negative (inward), effectively shutting down hyperpolarization and 

restoring the membrane to its resting state at the same time. 
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Figure 3: A circuit representation of model ([8]) of which a voltage-gating current replaces the 
leakage current in the HH model. The activation conductance and the gating conductance are of 
the same mathematical form but exactly opposite in the direction of the membrane potential. 


In addition, the sodium-gating parallel combination IV-cm\e, I = f^S^)+fQiV), is typically 
A^-shaped. The left knee point of the A^-characteristics is near the gating criticality Eq + and the 
right knee point is near the sodium criticality i^Na — This is because the tail part of each curve 
beyond its critical point is exponentially flat. Furthermore, the membrane resting potential is 
the intersection of the potassium IV-curve I = and the reflection of the sodium-gating 

IV-cmve I = -(/n,(^) + /g(^)) (so that f^{E^) + + fciEm) = 0). It is to the right 

side of Eq and below i?Na, and is always stable, see FiglU Negative conductance has always been 
a theoretical conundrum t ifTdl [TTl [181 iHl HI). For our model it is a mechanistic consequence to 
channel activation (|5]and voltage-gating ®. 

Conductance Adaptation and Circuit Model. Although each current’s IV-curve is accessible as an 
equilibrium current for each clamped voltage, in transient the voltage driving conductance g{V) = 
cannot be realized instantaneously during action potential generation and propagation. 
There is a time delay or time-course adaptation. We propose the following conductance adaptation 
following the idea of Hodgkin and Huxley (ifTI): 

^ - p) (7) 

dt 

with r being a time constant. Thus, instead of the IV-curve the corresponding current’s time 
course is given by / [t] = gp{t) (V(t) — E) with p{t) determined by the adaptation equation above. 

Putting all these assumptions together we obtain the following mathematical model for the 
squid giant axon and excitable membranes in general: 

' CV = -[g^niy - Ek) + -g^MV - ^Na) + -gyiv - Eq) - ly 
n' = — n) 

m' = — m) ^ ^ 

. h'= - h) 
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where Jgxt denotes any intracellular current other than the potassium, the sodium, the gating, and 
the capacitive currents. As for the conductance adaptation time constants, we assume the sodium 
channel and the gating channel share the same constant because they are structurally bonded 
together. The dimensionless variables n,m,h are not the same as of the HH model but are bor¬ 
rowed here to pay tribute to Hodgkin and Huxley’s seminal work. 

Best Fit of Model to Data for Parameter Estimation. A piece of mathematics only remains as a 
conceptual model for any physical process unless and until it is best-fitted to the process to fix a 
parameter point or a parameter range of the model. Otherwise different qualitative behaviors of 
the conceptual model would forever remain as unrealized and untested possibilities. To this end 
we re-sampled the action potential data of axon 17 from Fig. 12 of [HI and fitted our model to the 
data to see how well or badly the model performs. 

First we fix some parameter values and exclude them from the best-fit process because they are 
known. These are: the capacitance C, the two ions Nernst resting potentials Ek, £'Na, the intrinsic 
potassium conductance The membrane capacitance is C = l/iF per square centimeter. The 
Nernst potentials are extracted from [[T]| which had them shifted up by the Goldman-Hodgkin-Katz 
potential Er which was estimated by ifTTI to be around -47.5 mV, leading to E^ = —59.5 and 
^Na = 67.5 respectively. (As a result the HHAxonl7 data is shifted down by the supposedly 
membrane resting potential Er.) As for the intrinsic potassium conductance it is the value of 
g^^ (Ek)- As mentioned above it can be worked out to be g^^ = 0.0229 m.mho per square centimeter 
from [HI. 

The rest parameters are each given an estimated starting value and then best-fitted by a gradient 
search method. First we cannot do the same for the intrinsic sodium conductance as we did 
for the intrinsic potassium conductance g^ because we know now that Hodgkin and Huxley’s 
purportedly sodium JV-curve is an aggregate of the sodium and the gating channels. The estimate 
for the potassium conductance is reliable but not for the sodium conductance. Because of our 
estimation above on the order of magnitude for the intrinsic conductances of the potassium and 
sodium channels we will choose a starting value for the latter to be = 100. 

As for the gating parameter, we start the initial guess for the gating resting potential at Eq = 
—52.5mV, 5mV below the estimated membrane resting potential Er = —47.5mV. As for the in¬ 
trinsic gating conductance we will start it at g^. = 10, a value between g^ and 

As for the activation and gating range parameters we will start them at = 10, = 10, 5^ = 
8. This translates to the reversal of activation for the sodium channel at i^Na — b^^ = 49.5mV, a 
non-accessible and thus unimportant potassium criticality Ek — b^ = —69.5mV, and a gating 
‘threshold’ criticality Eq -f 5^ = —44.5mV. 

Last, for the adaptation time constants we use = 1/msec, and = 10/msec, making the 
sodium conductance adaptation one order faster than the potassium conductance. 

The HHAxonl7 data was obtained by Hodgkin and Huxley by an instantaneous depolariza- 
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Figure 4: (a) Following the initial guess from the text, the best-fitted parameter values for the 
eireuit model ([8]) are: -Ek = —59.5, = 0.0229, b-^ = 16.6, E'Na = 67.5, = 100, = 

18.4, Eg = -56, g^ = 9.3333, = 7.0667, C = I, = 0.8667, and = 10. Filled 

data points are used for the seeond-order error funetion. (b) The per-data-point error function 
in each searched parameter centered at the best-fitted point whose error value is 0.0082. (c) /17- 
curves of the best-fitted model, showing three equilibrium states with the left-most being the resting 
membrane potential Em- (d) A similar plot for the original fit of the HH model to the same data, 
and a best fit showing the voltage trace only (diamond marker) starting at the original parameter 
values of the HH model with searching parameters not shared with our model ([8]). The original fit 
error is 0.0595 and the best fit error is 0.0395, both are worse than the best fit of our model. 
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tion above the resting potential of the axon that translates to the initial voltage value V^(0) = 
—20.6707mV. The initial values for the rest of variables are n(0) = rn(0) = 

/i(0) = , whieh are fixed throughout the best fit proeess beeause they are good ap¬ 

proximations for all searehes at the initial time f = 0. 

A best fit of the model to the data is done by Newton’s line seareh method (e.f. [|20l llTTl f. the 
error funetion is defined to be 


N 

\\v{ti)-vn 


1/ 


where V is maximum of the absolute values of the voltage data, N is the number of data points 
ineluded for the error funetion. Thus £ measures the per-data-point relative error between the 
predicted value V{£) by the model and the observed value Vi at time fj. Fitted points include only 
these: the end points, the maximal and the minimal points, and the inflection points, constituting 
the so-called second-order fit. (The first order fit by definition would exclude the inflection points 
of the data.) The line search is to find the so-called best-fitted parameter values from the starting 
parameter values that gives rise to a local minimum point of the error function £. At any iteration of 
the search, the error is calculated at a discrete set of points from an interval of each parameter and 
the new starting parameter point is chosen if it defines the smallest error. The interval is centered 
at the parameter value with the radius of the absolute value of the parameter, i.e. either [0, 2p] or 
[2p, 0] depending on if the current parameter value p is positive or negative. FiglHshows both a best 
fit of the model and the original and a best fit of the HH model to the same data. This preliminary 
comparison suggests our model does better than the HH model. 

3. Discussion. Mathematical modeling is a process of falsification and refinement. The HH 
model has become a benchmark for the squid giant axon in particular and a template for excitable 
membranes in general. To replace its benchmark status by our circuit model, detailed comparisons 
between the two are needed and they are given below. 

Absolute Refraction and Anode Break Excitation. The HH model was hugely successful. Its suc¬ 
cess was supported among other things by its matching up two uncanny properties of Hodgkin and 
Huxley’s experiment data, the phenomenon of absolute refraction and the anode break excitation 
oscillation. The former occurs when a sudden initial depolarizing voltage is applied the membrane 
voltage decreases first before increasing or depolarizing as shown by both the axon 17 data and the 
models from Fig® This phenomenon was later identified to be the phenomenon of voltage-gating 
([Hill). It is validated by our model ([8]) as demonstrated in FigjSta). It shows not only model ([8]) ex¬ 
hibits the same property but also its affect by the intrinsic gating conductance as it is supposed 
to be. In contrast, the HH model exhibits the same but does not explain it. 

As for the anode break excitation, it refers to the phenomenon that after hyperpolarization and 
during the phase of re-polarization, the membrane voltage oscillates towards its resting potential 
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(a) (b) 

Figure 5: (a) The phenomenon of the so-ealled absolute refraetion by [[Tl for different depolarizing 
initial voltage with the same parameter values from FiglU for model dS]). Dotted eurves are the 
same plot exeept for a larger intrinsie eonduetanee = 15 for gating, (b) Anode break exeitation 
oseillations for small 0 < a < 1 for equations dH) with the same parameter values as (a) with 
= 15 exeept for the new additions a and = 0.5. 

rather than does so monotonieally. The HH model is eapable of this phenomenon but does not 
mateh up well in the time seale (Fig.22 of |[I1). To eapture this property with the eorreet time seale, 
we make a small modifieation to model dS]). We assume that a proportion of the gating eurrent 
{V — Eq) is subjeet to eonduetanee adaptation and the remaining proportion 
is subjeet to eurrent adaptation modeled by the last equation of the system below 

^ CV = -[gMV - Ek) + g^MV - ^Na) + ag^h{V - Eg) + (1 - aj/c - /ext] 

"E = — n) 

< m' = (9) 

E = _ h) 

. Ig =r^{goeE^-^^^'^<V-EG)-lG) 

where a : (1 — a) is the proportionality split for the two types of adaptation. Figl5](a) shows that 
if the proportion for eurrent adaptation is substantial (small a) but slow (small r^), anode break 
exeitation oeeurs with the same time seale as shown by the experimental data of |[I1 (Fig.22). We 
should note that the phenomena of absolute refraetion and anode break exeitation are notorious 
for other models to replieate, the FitzHugh-Nagumo model (UHl) or the Hindmarsh-Rose model 
( dull ) are two sueh examples. 

All-or-Nothing Action Potentials. Although both the HH model and our model dS]) mateh up 
Hodgkin and Huxley’s experimental data in many aspeets, many fundamental differenees remain. 
Figl^ shows some. The first eoneerns the all-or-nothing generation of aetion potentials. As it is 
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Figure 6: (a) The best-fitted aetion potential of ([8]) is projeeted to the i//K-spaee together with 
its various IV-curves, (b) Its eonduetanees as funetions of the time for the aetion potential, (e) 
and (d) Similar plot as (a) and (b) for the HH model. The dash-dot eurve is the refleetion of the 
sodium-leakage IV-curve I = —{fj^^V) + /l( 1^)) in (e). 

shown in Figlb^a), in addition to the membrane resting potential E'm, our model has two more 
equilibrium points: one to the right of Ey^, whieh is a saddle point, and one to its far right, whieh 
is typieally a unstable spiral. It is known from the theory of dynamieal systems that the resting 
membrane potential is always stable and the triple-equilibria configuration creates a saddle-node 
bifurcation point for the threshold of firing, giving rise to an all-or-nothing firing mechanism for 
action potentials. If the initial voltage is to the right of the middle saddle-node, an action poten¬ 
tial ensues, with a magnitude stretching passing the spiral focus at the least. In contrast, Figjbtc) 
shows that the HH model has only one equilibrium point which is the resting membrane potential 
Eya. The corresponding firing mechanism is by the way of a Hopf-bifurcation point. In theory, the 
action potentials are graded, not the all-or-nothing type as is supposed to be. That is, depending 
on how close the initial voltage is to the Hopf-point (the local minimum point of the reflection of 
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Figure 7: (a) With the same parameter values as from FigHa) for model ([8]) but foreed by an 
exeitatory square pulse I ext just below the aetion potential threshold whieh is about the maximal 
value of the gating eurrent. (b) The same plot exeept for an above-threshold square pulse, (e) A 
eyelie firing without the repolarizing resting potential if voltage-gating is absent (g^ = 0). (d) The 
sodium-gating IV-curves with varying gating range b^. Too small or too big a is bad for aetion 
potential generation. 


the joint sodium-leakage IV-cm\e I = -\- fj^iV))), the magnitude of an aetion potential 

ean vary from nothing to full, a phenomenon of the so-ealled eanard explosion ( Il^l24l f. 

There are three more notieeable differenees about the firing meehanism. For our model, the 
stability of the resting potential is due to the voltage-gating IV-curve, whereas that for the HH 
model is eritieally dependent of the presenee of the leakage eurrent Jl. Seeond, the re-polarization 
of the membrane for our model is due to the voltage-gating eharaeteristies as the sodium ehannel 
is elosing down when eharged moleeules return to the sodium gate pores. In eontrast, the re¬ 
polarization of the membrane by the HH model is earned out by the leakage ehannel. Third, our 
model prediets a persisting leaking eurrent when the membrane is at the resting potential as 


14 














































0 1 2 3 4 5 6 

Time (msec) 

(a) 



(b) 


Figure 8: (a) With the same initial eonditions, the same parameter values, and the same exeitation 
foreing as in FigjTtb) but with varying sodium aetivation range parameter (b) The eorrespond- 
ing IV-curve eonfigurations. Too small or too big a is not desirable for aetion potentials. 


shown in the inset of Figj^a) of whieh the primary eomponents are the gating and the sodium 
eurrents, eaneelling out eaeh other at the resting equilibrium. That is, the leakage is a eonsequenee 
of the voltage-gating rather than a eause by the HH model shown in the inset of Figj^e). 

As for the aetion potential at large, our model shows it traeks elosely around the sodium-gating 
eurve I = + f^iV)). Also the reversal of the potassium eurrent (trajeetory entering 

the left side of the potassium IV-curve) happens rather promptly after the potassium eurrent has 
peaked. In eontrast, the aetion potential of the HH model eomes elose only to two points of the IV- 
eurves, the Nernst potentials Ek, and with the rest playing little role to guide the trajeetory. 
And its potassium eurrent reversal is rather delayed after the potassium eurrent has peaked and the 
membrane has started depolarization. 

Voltage-Gating. One major differenee lies in the inelusion of voltage-gating into our model. As 
shown in Figj^b), voltage-gating is aetive during the initialization phase and the termination phase 
of an aetion potential. Both are for the purpose of preventing (aeeidental) firing of aetion potentials, 
and yet onee the firing threshold is exeeeded, the gating eonduetanee reeedes, permitting the mem¬ 
brane to depolarize. In eontrast, little ean be dedueed from the leakage eonduetanee time eourse 
(whieh is a eonstant of time) from the HH model as shown in Figj^d). FigjV] also shows more 
elearly what happens with the injeetion of an exeitatory intraeellular square pulse /ext(f)- Voltage¬ 
gating suppresses an aetion potential if the exeitation fails to elear the threshold, approximately the 
maximum of the gating eurrent. Otherwise an aetion potential is generated if the exeitation elears 
the threshold. The figure also shows that taking away gating {g^ = 0), the membrane resting po¬ 
tential is gone for losing its stability and the aetion potential is replaeed by a perpetual limit 
eyele. It also shows that if the gating range is too elose to the gating resting potential Eq (very 
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small b^) it is effectively the same as if there is no gating (g^ = 0). Also if the range of gating is far 
away (large b^) the threshold for action potential can become too high for intracellular excitation 


to clear. There is a ‘Goldilocks Zone’ from the resting gating potential for the range of gating. 


Ion Current Activation. FigjS] shows the effect of the sodium current activation range parameter 
b^^ on the generation of action potentials. If the sodium activation is too close to its Nersnt resting 
potential (small b^J it requires a very long excitation time for action potentials to arise. If the 
activation is far away (large b^J, it renders the voltage-gating useless. It also requires a ‘Goldilocks 
Zone’ from the resting sodium potential for the range of sodium activation. 

Model Overfit. Another major difference lies in the arbitrariness of the functional forms for the 
voltage-dependent conductances of the HH model. For example, for the function of the HH 
model, four parameters were actually required and fitted to have its final form obtained in |[ll 
because it is of this general form 


a2(exp(a3 — 04!/) — 1) 


Similarly, the function /3„ was fitted with two free parameters. Altogether, the HH model had 
13 more parameters for its fitting than our model (of which only the 6-parameters and the r- 
parameters are not shared with the HH model). Because the functional forms of the HH model and 
their parameters are rather arbitrary, its fit to the data can be construed as ‘overfit’ as with the case 
where arbitrary polynomials can fit but not explain any data. 

Traveling Action Potential. The term action potential in its original definition is referred to the 
uniform profile of the membrane voltage when an electrical pulse propagates from one end of the 
axon to another. The partial differential equations for such propagating pulse is derived by adding 
the axial diffusive current DV^x to the total current at each location to the patch model ([8]). The 
derivation follows the same treatment as in [HI. The spatially continuous model is 



( 10 ) 


with D being the axial diffusion coefficient of the axon. One can also derive a discrete version of 
the continuum model as follows 


CVl = -[g^ni{Vi - Ek) -f g^mfiVi - -f g^hfiVi - Eq) 

-Jext + - I",)] 



( 11 ) 


16 







50 

> 

0 

> 

-50 
0 

(a) (b) 

Figure 9: With the same parameter values as from FigHfa) for model (fTTl) . and d = 0.6 m.mho 
per square eentimeter, (a) shows the traveling pulse if the first node is injeeted by an intracellular 
current Jext = 35 yuA for the first 0.8 msec. Plot (b) shows the same except that the excitatory 
current is injected at the fourth node. 

with i = 1,..., n, di = dn+i = 0, and di = d, a constant for i = 2,..., n. The latter can also be 
used as a model for myelinated axon with i denoting the nodes of Ranvier ordered from one end of 
the axon to another. Figj^ shows the phenomenon of saltatory propagation by the discrete model. 

Spike Burst. We end this section with the phenomenon of spike burst for which the 

membrane potential of an excitable cell sustains a number of rapid spikes before hyperpolarizing 
and then re-polarizing near the resting potential of the membrane. Early models for bursting spikes 
include ll29l [30l [SB [32l [221. We will not fit our model to any specific data per se but rather 
highlight the simple mechanisms and configurations for spike bursting. The basic template for 
spike burst requires the all-or-nothing firing mechanism shown in Figj^a). The key difference 
between the one-spike action potential and the many-spikes burst is the time adaptation constant 
for the potassium conductance . For the former it has a slower time scale as shown in the figure 
where the sodium-gating reflection curve I = — (/^a (^) + /g (^)) strongly pulls the orbit to the left 
side of the potassium JE-curve / = If we speed up the potassium adaptation by increasing 

its time constant, the orbit will come down faster towards the potassium JE-curve / = and 

as a result oscillate around the spiral focus equilibrium point at the far right. That is, the middle 
saddle-node equilibrium can act to direct the orbit either to its left if the time constant is small 
or to its right to form spike burst if is large. Together, this forms the basic ingredients for spike 
burst. 

But in order to sustain multiple bursts, another ion current is needed to drive the bursting- 
spiking orbit between the two regimes, being separated and redirected by the middle saddle-node 
point. The dynamics is into the spiking regime when an orbit is directed to the right of the point 
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Figure 10: (a) Spike bursts of model ([8]) with the parameter values: £'k = —62, = 0.015, 

= 16, = 65, = 70, = 20, Eq = -50, g^ = 15, = 10, C = 1, = 3, 

Tj^^g = 80’ -^ci = —10, Vci = 10’ ^ci = bO, and = 0.02. The gating eurrent and the ehloride 
eurrent are magnified four times for a better view, (b) Exaetly the same parameter values for the 
3-D model (fT3l) . The IT-nulleline shown is its interseetion with Ici = 0. Dark dash line is the 
projeeted Jci-nulleline with Jci = 0 and the light dash lines mark the nulleline’s effeetive range for 
the spike burst. 


and into the quieseent regime when the orbit is direeted to the left of the point. For illustration 
purpose, we will inelude the ehloride ion Cl” eurrent to the membrane model. We will assume it 
has a negative Nernst potential Eqi < 0 above the resting membrane potential E^^, an exponential 
aetivation eonduetanee just like the other two ion speeies with moderate intrinsie eonduetanee ^gi 
and a eomparable aetivation voltage parameter b^^. Last, we will assume that it is eurrent-adapted 
and slowly so 0 < Tp; 1. The resultant system is as follows 

^ CV' = -[gMV - Ek) + g^MV - En.) + 9oKV - Eg) + /ci - /ext] 
n' = — n) 

< m ' = Tj^3,g ~ 

E = _ h) 

^ Jc/ = Tg, (V - Ecx) - lex) 

We also note that for many types of neuron and exeitable membranes their spike-burst dynamies 
are determined by sodium, potassium, and ealeium ions (Ca++) with the role of the sodium in the 
model above being replaeed by the ealeium (whose Nernst potential is usually higher still than that 
of the sodium ion) and the role of the ehloride replaeed by the sodium so that the region of spikes is 
always above the zero membrane potential V = 0. For mathematieal meehanisms whieh generate 
spike bursts, the prototypieal model above is nonetheless rather typieal and illuminating. More 
speeifieally, Fig JTCT a) shows a simulation of sueh a spike burst. One ean see that gating again plays 
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an important role for this behaviour, shutting down the burst when it is aetivated and permitting 
spiking when it reeedes. Also, it shows that the addition of the slow ehloride eurrent inereases 
during the spiking phase until sending the dynamies into the quieseent phase and deereases during 
the silent phase until sending the dynamies into the aetive phase. 

Spike burst is a 3-dimensional phenomenon just like aetion potential is basieally a 2-dimensional 
one. To demonstrate this point, let assume the sodium and gating dynamies is the fastest with 
~ oo so that their eonduetances equilibriumize instantaneously with m = and 

h = This eliminates two differential equations from (fT^ above to obtain the 3- 

dimensional system below 

^ CV = -[qMV - Ek) + _ E^,)+ 

g^e-^y-Ec)/^o iV - Eo) + Id - /ext] 

— n) 

. /ci' = - Ed) - Id) 

FigGSb) is the phase space of this system at the cross-section Id = 0. The l/-nullcline and 
the n-nullcline clearly show the all-or-noting firing configuration, having three intercept points. 
It also shows the spikes lie to the right side of the middle saddle-node of the 1/n-subsystem. As 
most part of the spikes lies to the right side of Jci-nullcline in which Ed > 0, Jci increases to 
become more positive. The outward current inhibits spiking, lowering the l/-nullcline until the 
orbit is caught to the left side of the middle saddle-node point and the spikes are switched off. 
But the quiescent phase lies to the left side of the /ci -nullcline in which Ed < 0, Jci decreases 
to become more negative, changing it into an excitatory current to drive the membrane into its 
spiking regime again. This is the basic 3-dimensional blueprint for continuous bursts of spikes 
in all higher dimensions. For exactly the same parameter values, the lower dimensional system 
(fT3]) and the higher dimensional counterpart (fT^ are all comparable except for possibly a different 
number of spikes. That is, if we project the higher dimensional bursts of spikes from Fig fTCT al to 
its 1/n-space, it should look like the phase portrait Fig fT^ bf of the lower dimensional counterpart. 

4. Concluding Remark. Our mathematical modeling of the squid axon began with a model for 
the sodium-potassium ion exchanger pump which led to a justification of the Nernst potentials for 
both ion species and to the modeling of the ion currents by battery-driven conductors. Following 
upon the voltage-dependent conductance finding of Hodgkin and Huxley we introduced the ex¬ 
ponential activation model for both opening and closing of ion channels. This activation model 
unifies two seemingly different types of activation in potassium and sodium channels. We also 
introduced the exponential voltage-gating model for the sodium channel which together with the 
sodium activation model automatically led to A^-shaped sodium-gating characteristics, solving the 
negative conductivity problem which has puzzled researchers of many generations. By incorpo¬ 
rating Hodgkin and Huxley’s time adaptation for channel conductances, our model is capable of 
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generating the all-or-nothing action potentials as a must-be consequence to voltage-gating as well 
as spike bursts with the all-or-nothing firing configuration. These voltage-gating related properties 
have never been demonstrated in any existing model. 

More specifically, our model enables the following narrative on action potential generation in 
mechanistic details unobtainable from the HH model. In the absence of external excitation the 
membrane is kept at rest by voltage-gating. When depolarizing excitation exceeds a threshold 
above the maximal voltage-gated current, the voltage-gating current drops exponentially to zero, 
allowing the membrane to depolarize which in turn opens up the sodium gates. The increased 
inward sodium current further depolarizes the membrane {dV/dt > 0) which in turn opens up 
not only more sodium gates but also potassium gates as well in the manner that their channel 
conductances grow exponentially with the depolarizing voltage. However, the inward sodium 
current must slow down as it is buffered by the potential difference between the membrane potential 
and its positive Nemst potential (/Na = fi'Na(^)(^ “ but in contrast the outward potassium 

current runs away exponentially from its negative Nemst potential. The outward potassium current 
must eventually catch up to the inward sodium current in magnitude so that the net ion and gating 
current becomes outward and the direction of the across-membrane voltage is reversed {dV/dt < 
0). This reversal activates the closing of both ion gates, closing up them exponentially fast. This 
in turn speeds up the downfall of the membrane potential. The hyperpolarizing potential may not 
pass its resting potential if the outward potassium current closes itself too quickly, giving rising 
to continuous spiking of some sort. Otherwise, the across-membrane potential must overshoot the 
voltage-gating potential {V < Eq), in which case the already-activated voltage-gating conductance 
is near its intrinsic conductance. From this point on the system must converge to the membrane 
resting potential whether or not the voltage overshoots it. This is because of two interplays of 
the currents. One, the voltage-gating current becomes negative below its reversal potential. Two, 
the sodium and potassium currents become so much smaller than the small inward voltage-gating 
current. So the latter becomes dominating, effectively halting membrane hyperpolarization and 
then reversing the direction of the hyperpolarizing membrane potential {dV/dt > 0), and bringing 
the membrane to its re-polarized resting state yet again. 

One conclusion seems obvious that our model ([8]) should replace the HH model as the basic 
template for excitable membranes for future researches in theoretical neuroscience. Another take¬ 
away from our result is that however complicated molecular biological processes may be it is not 
only possible but also imperative to model them mechanistically. Phenomenological models can 
fit but only mechanistic models can both fit and explain. Mathematical modeling is to find the 
mechanistic equation to which nature fits as a solution, as Newton demonstrated centuries ago. 
Richard Feynman believed that if we cannot lecture a theory to sophomore physics students then 
we must not understand the theory. Our model should pass his test. 
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